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We present the transverse momentum spectra and rapidity distributions of ir~ and Kg in Ar+KCl reactions 
at a beam kinetic energy of 1.756 A GeV measured with the spectrometer HADES. The reconstructed Kg 
sample is characterized by good event statistics for a wide range in momentum and rapidity. We compare the 
experimental tt~ and K% distributions to predictions by the IQMD model. The model calculations show that Kg 
at low tranverse momenta constitute a particularly well suited tool to investigate the kaon in-medium potential. 
Our Kg data suggest a strong repulsive in-medium K° potential of about 40 MeV strength. 
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I. INTRODUCTION 

Heavy ion collisions at relativistic energies in the 
SIS energy regime (E =1-2 A GeV) allow to cre- 
ate rather dense nuclear systems up to few times the 
saturation density and this provides a favorable envi- 
ronment for the study of in-medium hadron proper- 
ties. Within this context, expected medium effects on 
strange particles have been in the focus of nuclear re- 
action studies at SIS energies for the past two decades. 
The predicted appearance of a kaon condensate in 
compressed nuclear matter QjJ] with its consequences 
for the understanding of neutron star evolution ||2|] has 
emphasized the quest for a quantitative determination 
of the Kaon-Nucleon/Nucleus potential. So far, partic- 
ular efforts have been devoted to the production of K + 
and K~ mesons and have provided phase space distri- 
butions, integral yields and angular distributions for a 
wide range of energies and collision systems. 
The systematics of the experimental K + observables 
and in particular the results from sideward JH] and out- 
of-plane [J] flow analyses suggest a repulsive kaon- 
nucleus potential. Data from proton-induced reactions 
support a moderately repulsive potential for K + of the 
order of 20 MeVfl 0] in agreement with theoretical 
calculations |p7l— TlQtl . 

Results concerning the K~ lUt] are hampered by the 
low statistics available, and no settled quantitative con- 
clusions could be drawn so far about the strength of the 
attractive potential. 

Neutral kaons can shed additional light on the under- 
lying question. They have the advantage that possible 
medium effects are not obscured by the Coulomb inter- 
action. Similar to K + , the K° in-medium potential is 
expected to be repulsive at these energies as can be in- 
ferred from recent results extracted from pion-induced 
reactions [11]. Indeed, the comparison of the K° s mo- 
mentum distribution in the tt~ + C and tt~ + Pb reac- 
tion points to the existence of a repulsive KN potential 
of 20 ± 5 MeV at a normal nuclear density. 
In this work we report on results for tt~ and Kg ex- 
tracted from Ar+KCl reactions at 1.756 A GeV. The 
high statistics data sample covers almost the entire 
phase space and allows a detailed analysis of the low- 
momentum component. For the first time, p t distribu- 
tions for Kg down to 50 MeV/c for the whole rapid- 
ity range could be measured. The study of the low- 
momentum region is well suited to access the K° po- 
tential in the nuclear medium since there repulsive ef- 
fects are expected to show up in a more pronounced 
way. Moreover, the spectral shape of the p t distribu- 
tion allows quantitative conclusions, without requiring 
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any absolute normalization necessary for descriptions 
which rely on measured yields only. The obtained K% 
data are compared to our results obtained for K + B12I1 
and to theoretical calculations by the IQMD (Isospin 
Quantum Molecular Dynamics) transport model lfl3ll 
leading to an estimate of the strength of the repulsive 
K° in-medium potential. 

Our paper is organized as follows. Section II sum- 
marizes the experiment and section III addresses as- 
pects of the particle identification method. This sec- 
tion presents also transverse mass spectra and rapidity 
distributions for ir~ mesons including a comparison to 
IQMD results. Section IV is devoted to the Kg recon- 
struction procedure based on ir~7r + pair decays. In 
section V we contrast the obtained K s distributions to 
K + data measured in the same reaction and compare 
them to results of IQMD simulations. In this section 
we discuss the findings concerning the K° in-medium 
potential. We close with a summary in section VI. 



II. THE EXPERIMENT 



The experiment was performed with the High 
Acceptance Di-Electron Spectrometer (HADES) at the 
heavy-ion synchrotron SIS at GSI Helmholtzzentrum 
filr Schwerionenforschung in Darmstadt, Germany. A 
detailed description of the spectrometer is presented in 

H. 

HADES consists of a 6-coil toroidal magnet centered 
on the beam axis and six identical detection sections 
located between the coils and covering polar angles 
from 18° to 85°. In the measurement presented here, 
the six sectors comprised a gaseous Ring-Imaging 
Cherenkov (RICH) detector, four planes of Multi-wire 
Drift Chambers (MDCs) for track reconstruction and 
two Time-of-Flight walls (TOF and TOFino), sup- 
plemented at forward polar angles with Pre-Shower 
chambers. For each sector, the TOF and TOFino/Pre- 
Shower detectors are combined to a Multiplicity and 
Electron Trigger Array (META). 
A fg+Ar beam of ~ 10 6 particles/s was incident on 
a four-fold segmented KC1 target with a total thick- 
ness corresponding to 3.3% interaction length. A fast 
diamond start detector located upstream of the target 
was used to determine the interaction time. The data 
readout was started by a first-level trigger (LVL1) de- 
cision, requiring an observed charged-particle multi- 
plicity MUL > 16 in the TOF/TOFino detectors, 
accepting approximately 35% of the nuclear reaction 
cross section. This centrality selection translates into 
an average participants number A part = 3 8. 5 ±2. 7, as 
described in IHflBll . 
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III. PION IDENTIFICATION 



The reconstruction of K° s —> n~ir + decays requires 
a clean and unambiguous pion identification within 
the high multiplicity track ensemble of each event. 
Charged particles fire wires in the MDCs in front and 
behind the magnetic field and one or more hits in the 
META detector. In the so-called cluster finder soft- 
ware Il4 UM track-segments are formed using the in- 
formation from the two pairs of MDC planes. In the 
track-segment fitting procedures possible trajectories 
through the two track segments and the META hits 
are calculated. For a detailed description of the track- 
ing procedure see lfl2l [nil . After proper correlation 
of track segments and META hit points, particle mo- 
menta were calculated with a Runge Kutta integration 
of the trajectory in the magnetic field. The graphi- 
cal cuts utilized to select the tt^ 1 sample on the base 
of the dE/dx versus momentum distributions are the 
same as shown in 11211 . This selection allows to extract 
a high purity pion sample. Muons from pion decays 
(tt ± -)• p.* + v^, BR = 99, 9%, ex = 7.8 m) inside 
the HADES spectrometer constitute only a small frac- 
tion of the tracks and are mostly misidentified as pions 
due to their small mass difference. Hence, the effective 
losses are in the order of few percent only, as shown in 

a. 

The efficiency of the MDC - dE/dx cut and the purity 
of the selected sample have been extracted from ex- 
perimental data, selecting particles fully reconstructed 
in the spectrometer and identified by the META as a 
reference and then verified using simulations, apply- 
ing to the simulated tracks the same selection criteria 
as for the real data. The average efficiencies, defined 
as the fraction of real pions surviving all the cuts are 
~ 90%. The purity of this sample also corresponds 
to ~ 90%. While these values are rather independent 
of particle momentum for ir~ mesons, a significant re- 
duction of the purity is observed for tt + with momenta 
p> 400MeV/c due to contamination by misidentified 
protons. However, the Kg reconstruction is only very 
weakly affected by the purity of the high momentum 
7r + sample. The acceptance of the HADES spectrome- 
ter and the efficiency of the tt^ reconstruction were de- 
termined by a full-scale GEANT3 [19] simulation. A 
sample of pions generated with a flat distribution in ra- 
pidity (—0.75 < y c .m. < 0.75) and transverse mass (0 
< m t — m v < 550 MeV/c 2 ) was propagated through 
the whole spectrometer yielding the geometrical ac- 
ceptance from registered detector hits. The detection 
efficiency of the pion sample was obtained from the 
ratio of all emitted tracks inside the acceptance to fully 
reconstructed tracks subjected to the same selection 
criteria as for the experimental data. 
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FIG. 1: (Color online). Reduced transverse mass distribution 
for different rapidity bins for the backward hemisphere for 
7r~ together with predictions by the IQMD model (dashed 
lines) and the fit (full lines) of the experimental data (sym- 
bols). 



IV. PION RESULTS 

Reduced transverse mass distributions for tt~ are 
shown in Fig. Q] for different center of mass rapidi- 
ties (y c . m .=y-y(c.m.), y(c.m.)= 0.86). The data are 
corrected for the efficiency and the spectrometer ac- 
ceptance and are normalized to the number of LVL1 
triggers. The representation per transverse mass and 
rapidity unit divided by is chosen to ease the com- 
parison with a Boltzmann distribution. A fit according 
to 



1 d 2 M 
to j dm t dy 



C{y) exp 



(mt - 777q)c s 

T B (y) 



(1) 



has been applied to the m t — distribution in two 
adjacent m t intervals. The pions mainly stem from the 
decay of the A(1232) resonance and are produced at 
different stages of the collision with varying 'hardness' 
as time proceeds. Hence, the measured it spectrum is 
a superposition of different A generations. This ex- 
plains why a single Boltzmann fit (Q]) is not sufficient 
to match the data. It was found that a simultaneous 
two-slope Boltzmann fit with the transverse mass in- 
tervals chosen as < m t — < 180 MeV/c 2 and 
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FIG. 2: (Color online). (Panel a) Rapidity dependence of 
slope parameters obtained by fitting the n~ m t distributions 
from experimental data (full symbols) and IQMD simula- 
tions (dashed curve) in the range < m t — < 180 
MeV/c 2 . (Panel b) Slope parameters obtained fitting the 
same distribution in the range 180 < m t — m-n < 500 
MeV/c 2 . The values of the effective slopes T e gi/2 have been 
obtained employing function (O. 



180 < mt — < 500 MeV/c 2 , reproduces the data 
adequately. The two m t ranges have been first fitted 
with two independent Boltzmann distributions yield- 
ing the start parameter values for the simultaneous fit. 
For the latter, the start values were allowed to vary up 
to a maximum of 10%. The resulting Boltzmann dis- 
tributions are plotted in Fig.[T]as solid lines. 
As already shown in |20], measured pion spectra could 
be reproduced reasonably well by IQMD calculations 
Il2ll427ll . The present HADES data are compared to 
IQMD calculations assuming a hard cut on the impact 
parameter b < 6 fm which corresponds to a num- 
ber of participants of 38.5±4.6. The ir~ m t distri- 
butions obtained for the experimental data (symbols) 
are shown in Fig. Q] for different y c . m . bins together 
with the IQMD calculations (dashed lines) assuming 
the centrality selection b < 6 fm. The same fitting pro- 
cedure has been applied to the IQMD data points. The 
obtained inverse slope parameters T± t 2 for experimen- 
tal and simulated data are shown in Fig. [2] as a function 
of rapidity y c m .. We observe a systematic difference 



of 10-15% between the slope factor T\ determined in 
the forward and backward rapidity hemispheres which 
is not present in T%. The agreement between data and 
model is of similar quality. 

The full curves shown in Fig. [2] represent a fit to the 
experimental data according to 



Tb(v) 



T, 



off 



COSh(y - y c .m.) 



(2) 



that allows to extract the effective slopes for the pion 
in the two m t ranges. The obtained values T ffi 2 are 
shown in the two panels of Fig. [2] 
Transverse momentum spectra of pions in Ar+KCl col- 
lisions at a very similar energy (1.808 A GeV) have 
been published as kinetic energy spectra at 90° in the 
centre-of-mass frame ll28ll . which corresponds to mid- 
rapidity. The comparison with our data after proper 
normalization to the same number of participant nu- 
cleons shows agreement within the quoted errors. 
The functions resulting from the fits to the pion 
transverse-mass distributions in the different rapidity 
bins in Fig. Q] are integrated over the interval < 
m t — < oo in order to obtain the yields as a func- 
tion of rapidity. This integration method is justified 
by the large coverage of the experimental m t — 
distribution and the good agreement between the data 
and the fit function. The resulting rapidity distribution 
of 7r~ is displayed in Fig. [3] together with the values 
obtained integrating the IQMD distributions. The two 
distributions are in agreement within 15%; the largest 
deviation shows up in the mid-rapidity range. The to- 
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FIG. 3: (Color online). Rapidity distribution of ty~ exper- 
imental data and reflected data points around mid-rapidity 
(full and empty triangles) together with the result from the 
IQMD calculation (dashed curve). The full curve shows a 
Gaussian fit of the experimental data. 



tal yield of negatively charged pions was calculated 
from a Gaussian fit to the dN/dy distribution shown 
in Fig. [3] The systematic error on the pion multiplic- 
ity has been estimated varying the boundaries used as 
start values for the transverse mass fits in the range 



from 120 to 300 MeV/c 2 and repeating for each step 
the same fitting procedure described above. We obtain 
M exp (7T-) =3.9±0.1±0.1andMi Q MD(>-) = 3.7 
which are compatible within the error. Any system- 
atic error derived from the comparison of the backward 
and forward distributions was found to be smaller than 
the statistical error. These results have been confirmed 
by an independent analysis in 01511 . Taking as refer- 
ence the pion multiplicities estimated in Il20ll for the 
Ar+KCl system at 1.756 GeV, agreement is achieved 
as well. 

A comparison with earlier data and the IQMD model 
which is independent of the trigger details and im- 
pact parameter selection is possible by normalizing 
the pion yield to the number of participating nucleons. 
For the data presented here we find (N^)/A p&r t = 
0.101 ±0.003±0.008 (where the last error contains the 
uncertainty of the calculation of the A par t value with 
UrQMD, « 7%). From IQMD we get (N w )/A part = 
0.105 ± 0.015. From i§ for (N„) = 5.6 and A part 
=60 we obtain 0.093. 



V. K% RECONSTRUCTION 



The K s meson (mean decay length ct= 2.7 cm) de- 
cays into a 7r + + 7r~ pair with a branching ratio of 
69%. After selection of events with two identified 
pion tracks, the invariant mass was calculated for each 
7T + — 7r~ combination. In order to reduce the combina- 
torial background from uncorrected tt + — tt~ pairs and 
hence to enhance the purity of the Kg signal, various 
cuts on characteristic geometrical distances have been 
applied: (1) the minimum distance between the two 
pion tracks (d 7r +_ w - < 10 mm), (2) the distance of 
closest approach to the primary vertex for the two pion 
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FIG. 4: Invariant-mass distribution of 7r + — tt~ pairs (panel 
a). The combinatorial background (shaded area) is obtained 
by the mixed-event technique. The background- subtracted 
distribution (panel b) shows clearly the Kg meson signal 
(grey area with a Gaussian fit). 



The width of the reconstructed signal and the signal-to- 
background ratio depend on the rapidity bin and vary 
in the range 7MeV/c 2 < u K o < 12.4MeV/c 2 and 
0.3 < S/B < 2.0, respectively. A total of about 
65.700 K s mesons are identified in an interval of ±3o~ 



and (3) the distance between the primary reaction and around the fitted mass peak. 



secondary decay vertex (d(K s — V) > 30 mm). The 
method is discussed in detail in ref. 11711 . Fig.|4](top) 
shows the resulting invariant-mass distribution of all 
7T + — 7r~ pairs. The peak corresponding to the K s 
signal is clearly visible on the top of the background 
which is reproduced using the mixed-event technique 
(dashed area in the top of Fig. |4j. For the background 
analysis, only those events were chosen for which the 
individual pion tracks originate from the same target 
segment of the four-fold KC1 target stack. 
The resulting Kg signal after background subtraction 
is shown in Fig. [4] (bottom) and it can be fitted by the 
sum of two Gaussian distributions. A rather sharp dis- 
tribution that contains most of the yield and a broader 
distribution that corresponds to the case where at least 
one of the two pion has undergone multiple scatter- 
ing. The Gaussian fit in Fig.|4]shows the result for the 
sharper distribution and corresponds to the following 
mean value and dispersion for the reconstructed mass: 
(m K o) = 492.5 MeV/c 2 and (cr K o) = 9.3MeV/c 2 . 



The phase space distribution of the measured K° s is 
shown in Fig. [5] as a function of the reduced transverse 
mass in t — m K o and center of mass rapidity y c . m .. 
From the plot it is evident that the low transverse mo- 
mentum region is covered with significant statistics for 
the whole rapidity range. 

For the quantitative analysis, acceptance and ef- 
ficiency corrections have been applied to the recon- 
structed Kg signal in a similar way as described for the 
single w tracks in section [Til] A K s event generator 
with a flat distribution in rapidity (—0.75 < y c .m, < 
0.75) and transverse mass (0 < m t — m K o <900 
MeV/c 2 ) has been used as the input of a full-scale sim- 
ulation in order to evaluate the geometrical acceptance 
and the reconstruction efficiency. The average geomet- 
rical acceptance for K s amounts to 20-25% 1 1711 . the 
reconstruction efficiency to ~ 5 — 10%. The latter is 
shown in detail in Fig.|6]as a function of the transverse 
momentum p t for different rapidity bins. In addition, 
the branching ratio of the decay K s — >• tt + + n~ was 
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FIG. 5: (Color online). Distribution of reconstructed Kg 
yield as a function of the subtracted transverse mass m t — 
m K o and the center of mass rapidity y c . m .. The color code 
refers to the number of kaons per bin. 
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FIG. 6: K5 reconstruction efficiency as a function of the 
transverse momentum p t and for different rapidity bins. The 
geometrical acceptance is not taken into account. 



VI. k° s results 

A. Transverse Mass Spectra 

The Kg transverse mass spectra obtained after back- 
ground subtraction and correction for acceptance and 
efficiency are plotted in Fig. |7](open symbols) for var- 
ious y c m . bins together with the K + spectra (full sym- 
bols). Since in isospin symmetric heavy-ion reactions 
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FIG. 7: (Color online). Transverse mass spectrum for differ- 
ent center of mass rapidity bins, corresponding to the back- 
ward hemisphere only, for K s (empty symbols), K + multi- 
plied by a factor 0.5 (full symbols) and fits assuming a Boltz- 
mann parameterization (dashed lines). 



the yields of K + and K° should be the same and since 
the K~ yield is negligible compared to the K + yield we 
have K+ + K- = K° + K° = K° s +K a L = 2-K% For 
this reason, the data for the K + mesons (full symbols) 
in Fig. [7] have been multiplied by a factor 0.5. The 
quantitative comparison of these experimental spec- 
tra, that are found to be in good agreement, is a valu- 
able cross-check of the analysis. Neglecting final state 
Coulomb interactions, the dynamics of the interaction 
between K + and K° mesons and the nuclear medium 
should result in similar kinematic distributions. 
Using the Boltzmann parameterization ((T), an inverse 
slope parameter Tg(y) can be determined as a func- 
tion of the rapidity. The dashed lines in Fig. [7] shows 
the Boltzmann fit-functions. The resulting Tb{u) val- 
ues as a function of the rapidity obtained from the fit of 
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the K° data are shown in Fig. [8] together with the re- 
sults obtained for the K + [12] and by the IQMD simu- 
lations. The two curves in Fig. prefer to two different 
scenarios: 1) no in-medium potential (dotted lines in 
Fig.|9]l and 2) a repulsive potential of 40 MeV (dashed 
lines). Fitting the Ts(y) distribution with the function 
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FIG. 8: (Color online). Rapidity dependence of slope pa- 
rameters obtained by fitting the K5 and K + m t distribu- 
tions from experimental data (empty and full symbols re- 
spectively). The results by IQMD simulations employing a 
repulsive potential (dashed curve) and no potential (dotted 
curve) are shown as well. The errors are statistical only. 
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(|2j, the parameter T e g can be extracted. It represents 
the inverse slope at mid-rapidity and corresponds to an 
effective temperature at the kinetic freeze-out stage. 

The averaged T cS for is found to be 92 ± 2 MeV 
which agrees within the errors with the findings for K + 
in the same data set lfl2ll . 



FIG. 9: (Color online). Rapidity distribution of Kg (empty 
triangles) and K + mesons (full triangles from 1 12]) together 
with the distributions obtained with the IQMD calculation 
assuming a repulsive (dashed line) and no (dotted line) po- 
tential for two different normalization procedures (panel a 
and b, see text for details). The vertical boxes show the sys- 
tematic errors associated to the K + data. 



B. Rapidity Distribution 

The fitted K° s invariant transverse mass distribu- 
tions are integrated within the whole interval < 
m t — m K o < 00 in order to obtain the meson yield 
per rapidity unit. To evaluate the precision of this in- 
tegration, the integral of the fitting function has been 
compared to the integrated data points in the range 
< m t — m K o < 600 MeV/c 2 . We find an agreement 
between 0.5 and 2.5% for all rapidity bins, except for 
the interval —0.65 < y c . m . < —0.55 for which the dif- 
ference amounts 6.5%. 
The resulting rapidity distribution is shown in Fig. [9] 
(empty symbols), together with the respective K + 
yields divided by a factor two (full symbols). The 
systematic error of the K° s yield has been estimated 
by varying the geometrical cuts and the track quality 
selection in several combinations. The systematic er- 
rors shown for the K + distribution have been calcu- 
lated as described in [12]. Both experimental distribu- 
tions overlap within the error bars. Appropriate inte- 



gration yields a total Kg multiplicity of (1.15 ± 0.05 ± 
0.09) • 10~ 2 as compared to a total K + multiplicity 
of (2.8±0.2±0.1±0.1)- 10~ 2 [12], where the first and 
second (third) error are the statistical and systematic 
respectively. The differences in the tails of the rapidity 
distributions shown in Fig. [9] lead to slightly different 
47T yields. Nevertheless the multiplicities, extracted 
with two independent analyses, agree within the error, 
confirming the quality of the K s data. 
Together with experimental distributions the IQMD re- 
sults for the same two scenarios presented in Fig. |8]are 
shown. In the top panel of Fig. [9] the IQMD calcula- 
tions (b < 6 fm) without any normalization are shown, 
while in the bottom panel the two curves are shown 
after a normalization to the total area of the experi- 
mental distribution. The IQMD calculation used here 
corresponds to the standard setting with a — 1.0 (see 
section IVTCl already shown in lflq|29tl . 

Although the standard parametrization used in the 
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IQMD calculations above has been successfully used 
for comparison to KaoS and FOPI data and thus sup- 
ports the results of this comparison, a direct conclusion 
on the optical potential from absolute yields is prema- 
ture. It was shown in B3011 that transport models us- 
ing opposite assumptions on the optical potential but 
different parametrizations of poorly known production 
cross sections were able to reproduce the same rapid- 
ity distribution of kaons in central Ni+Ni events at 1.93 
AGeV measured by FOPI and KaoS. 

The calculations with and without potential have 
indeed very similar shapes (bottom panel of Fig. [9). 
The rapidity distribution contains the integrated 
information about pt and assuming that the effect of 
the potential is momentum dependent, the differential 
study of the p t distribution can add important informa- 
tion. The study of the p t distributions between 50 and 
800 MeV/c constitutes the key clue of our approach to 
determine the in-medium K° potential, as described 
in the following section. 



C. Comparison with IQMD 

As shown in section [TTTJ the IQMD calculations are 
in reasonably good agreement with the pion spectra 
measured for the Ar+KCl reaction at 1.756 A GeV 
as far as the slopes are concerned, but differ up to 
15% in the absolute yield. These findings impose a 
caveat on the absolute normalization of the IQMD sim- 
ulations on the experimental data. Starting with the 
same centrality selection used to obtain the pion spec- 
tra (b < 6/m), a set of calculations with the IQMD 
model 11311 has been carried out, employing a repul- 
sive K°-nucleus potential of varying strength. 

According to JU], K + = su is a "good" quasi- 
particle with narrow width; the same is expected for 
K° = sd. Its dispersion relation may be written as 
cu 2 = m 2 K + k 2 + II with cu(fc) as energy (mo- 
mentum) and II as real part of the self-energy. The 
latter one accounts for the influence of the ambient 
medium, thus leading to an effective in-medium mass 
to* = w(k = 0). In agreement with the low-density 
theorem, the K° effective mass should increase in nu- 
clear matter, as for the K + . Another way to state this 
is to say that the kaons interact with surrounding nu- 
cleons via a Schrodinger-type potential. An increase 
of m* is thereby related to a repulsive total poten- 
tial. The easiest parameterisation is the linear Ansatz 
to* = m p + U(a)p/po where p(po) is the nuclear 
matter (saturation) density. For the IQMD settings 
lfl3tl . the parameter a is related to the potential U via 
U(a) ~ U + U'a with U « 0.8 MeV and U' 38 
MeVH. 

In order to evaluate more quantitatively the compar- 
ison between IQMD and experimental data, we have 
first focused on the mid-rapidity p t distribution and 



x10 -6 -0.07<y e m <0.07 




200 400 600 800 .5 1 1.5 



p t [MeV/c] « 

FIG. 10: Color online. Left panels: p t distributions of Kg 
at mid-rapidity (full symbols) compared with different cal- 
culations by the IQMD model. The different solid curves 
correspond to a variation of the parameter a, which couples 
directly to the potential strength. The comparison is shown 
for three different normalization procedures (see text for de- 
tails). Right panels: Normalized x 2 distribution as a func- 
tion of the parameter a extracted from the comparison of the 
IQMD calculations to the experimental data for three differ- 
ent normalization procedures (see text for details). The insets 
show a zoom around the minimum of the distribution. 



compared systematically the experimental data with 
different IQMD calculations obtained varying the pa- 
rameter a. The left panels of Fig. [10] show the K s p t 
distribution at mid-rapidity, together with the results 
of the IQMD calculations for different values of the 
parameter a. The three panels correspond to three dif- 
ferent ways of normalization of the IQMD calculations 
to the experimental data. Panel (a) in Fig. [TOl shows 
the case where the IQMD curves for b < 6 fm have 
not been normalized to the data, panel (b) shows the 
case where all the IQMD curves have been normalized 
to the high-pt tail (pt > 475 MeV/c) of the experi- 
mental distribution and panel (c) shows the case where 
the area underlying each curve has been normalized to 
the integral of the experimental distribution. The error 
shown for the experimental data contains also the sys- 
tematic contribution. This contribution has been eval- 
uated varying the cuts for the Kg selection for 14 dif- 
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ferent combinations. A x 2 analysis has been carried 
out applying a best fit of the different IQMD curves 
to the experimental data. The result corresponding to 
the three normalization procedures are shown in pan- 
els (a'), (b') and (c') of Fig. [TOl where the error bands 
include both the statistical and systematic contribu- 
tions. The inlets in the three panels show a zoom on 
the minimum region. The minimum of the x 2 distribu- 
tion that corresponds to an optimal set of cuts for the 
Kg is obtained for the following a values respectively: 
1.13 - 0.12, 1.37 - 0.2 and 1.34 - 0.17. 
The asymmetric error on a has been evaluated taking 
the minimum x 2 value +1 for each of the 14 differ- 
ent x 2 distributions and reading the corresponding a 
value on the left side of the minimum. The maxi- 
mal deviation of this value from the minimum of the 
dashed curves has been assigned to the asymmetric er- 
ror. Since IQMD calculations corresponding to a val- 
ues higher than 1.5 lead to unphysical results, it is not 
possible to calculate the upper value of the a parame- 
ter. 

One can see that even considering the three normaliza- 
tion methods simultaneously the lowest limit for a is 
1, which correspond to a minimum repulsive potential 
of 38.7 Me V. 

It has also to be mentioned that the IQMD calculations 
quoted in [29], which interpret the whole K + KaoS 
systematics in favor of a soft EOS but do not draw any 
conclusions on the potential, used a — 1.0 as well. 
This rather high value of the extract ed p otential is 
larger that the results reported in S HL S [Of l22l , which 
reported a value of 20 ± 5 MeV. However, the potential 
values extracted so far from K s data are derived from 
7r~ +A reactions HXlfl - In this case the pion absorp- 
tion happens on the nucleus surface, so that many of 
the produced kaons do not travel through the nucleus. 
Furthermore, the HADES data deliver higher statistics 
and accuracy in the measurement of the low p t range, 
which is more sensitive to the potential effects. Also, 
the results extracted from K + data in proton-induced 
reactions 1 32] test subnormal nuclear density and can 
not be directly compared with the data from heavy- 
ion collisions where average densities of 1 .5-2 po are 
reached. In overall the nuclear environment resulting 
in the Ar+KCl reaction could lead to a stronger repul- 
sive potential for K°. 

On the other hand, as extensively discussed in fioll33ll . 
the different transport codes that have been employed 
to interpret the available data on kaons differ substan- 
tially in their implementation, in their elementary pro- 
duction cross-sections and in their treatment of the pa- 
rameters connected to the strength of the potential. To 
get a consistent picture, the HADES Kg data can be 
described by other models and the interpretation of the 
K + data previously measured in heavy-ion collisions 
yD might be revised using updated versions of theoret- 
ical models and elementary cross-sections ll33tl . 

Taking as a reference the IQMD simulation cor- 
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FIG. 11: Color online. Left panels: pt distribution of the ex- 
perimental Kg data (full triangles) together with the yields 
calculated by the IQMD model including a repulsive K°- 
nucleus potential of 46. 1 MeV (dashed curves) and without 
potential (dotted curves) for different rapidity bins. Right 
panels: ratio between the calculation by the IQMD model 
and the experimental data as a function of pt for different 
rapidity bins. 



responding to a=0 and a=\.2, we have compared 
them to the experimental p t distributions as shown in 
Fig-El The left panels of Fig.fTTIshow the K° s p t dis- 
tributions for different rapidity bins, together with the 
IQMD calculations assuming either no repulsive po- 
tential (dotted curves) or a potential of about 46. 1 MeV 
(a — 1.2, dashed curves). For this comparison the sim- 
ulations are normalized to the high-momentum tail of 
the experimental distribution starting from p t = 475 
MeV/c, since in this region the effect of the repulsive 
potential should be negligible. The displayed errors are 
only statistical. On the right panels of Fig. [TTJ the ratio 
of the simulated to experimental yield is displayed for 
the two cases with and without a repulsive potential. 
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The error band displayed for the ratios are again only 
statistical. 

It is very interesting to notice that the IQMD simula- 
tion assuming a repulsive potential agrees rather well 
with the experimental data for all the rapidity bins. The 
calculations without the potential, on the other hand, 
overestimate the experimental data, especially in the 
low transverse momentum regime (< 400 MeV/c), as 
is clearly visible in the ratio plots. This behavior is 
slightly more evident at mid-rapidity but persists with 
the same trend for the whole rapidity range. 
We note that the p t spectra exhibit the most sensitive 
dependence on the K°-nucleus potential. The rapidity 
distribution is sensitive to variations of the quantity a 
as well, but this effect is less pronounced and the con- 
clusions drawn on the potential strength are strongly 
dependent upon the normalization of the simulated to 
experimental data. Hence the accurate study of the p t 
spectra delivers a new more powerful tool to extract 
quantitative information. 



VII. SUMMARY 

Rapidity and transverse momentum spectra of tt~ 
and K° s mesons produced in Ar+KCl collisions at 
1.756 A GeV and measured with HADES have been 
investigated. For the first time in this energy regime, 
Kg spectra have been measured with high statistics and 
precision in almost the full phase space and, in partic- 
ular, down to low momenta (p m in ~ 50 MeV/c). 
The 7r~ data have been shown and compared quanti- 
tatively to IQMD calculations, using an absolute nor- 
malization based on the selection of the impact param- 
eter b < 6 fm in the simulations and a centrality se- 
lection on the experimental data corresponding to the 
most central 35% of the total cross section. The 7r~ ra- 
pidity density distribution shows a quantitative agree- 
ment with the IQMD model within 15%. 
For Kg mesons the transverse mass and the rapidity 
density distribution are compared to previously pub- 
lished K + data and agree well with them, both in shape 
and yield. The comparison of the Kg. rapidity distribu- 
tion with IQMD calculations does not seem to be solid 
enough to extract reliable information about the poten- 
tial, since the normalization is not certain and the cal- 
culations corresponding to different values of the po- 
tential deliver curves with the same shape. The Kg p t 
distributions have been found to be a better observable 
and have been compared to calculations by the IQMD 
model assuming different strengths of the K°-nuclear 
medium potential. This comparison supports the exis- 
tence of a rather strong repulsive potential of about 40 
MeV. These data are now available for further studies 
via the available transport models. 
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